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O ■ Abstract 

(N ' 

We study the behavior of eigenvalues of matrix pv + Qn where P/v and Q n 
C$ • are two N-by-N random orthogonal projections. We relate the joint eigen- 

value distribution of this matrix to the Jacobi matrix ensemble and establish 
the universal behavior of eigenvalues for large N. The limiting local behavior 
of eigenvalues is governed by the sine kernel in the bulk and by either the 
■ Bessel or the Airy kernel at the edge depending on parameters. We also study 

an exceptional case when the local behavior of eigenvalues of P/v + Qn is 
^ . not universal in the usual sense. 

<S . 

1. Introduction and Notation. 



q 
in 



Let A and B be two complex Hermitian matrices, perhaps random, and consider matrix 
X = A + UBU* where U is a uniformly-distributed unitary matrix. What can be said about 



> 

m 

Q\ • the distribution of eigenvalues of X if the size of the matrices is large? Especially, what can 



be said about the properties of this distribution on the local scale? 

It is known that if B is a random matrix with Gaussian entries, then the eigenvalues of A + 
B can be described as a determinantal process (Brezin-Hikami Q). Moreover, if in addition 
A is a Wigner matrix (that is, a Hermitian matrix with independent upper-diagonal entries) and 
if the entries of A and B are comparable in size, then the local properties of matrix A+UBU* 
are the same as if this matrix were from GUE, Gaussian Unitary Ensemble, (Johansson lfl3l ). 
On the other hand if B is a multiple of the identity matrix then the eigenvalues of A + UBU* 
are obviously simply the shifted eigenvalues of the matrix A. 

The main question at stake here is that of the universality of the eigenvalue distribution of 
A + UBU*. The universality hypothesis claims that for a large class of matrices the local 
correlations of eigenvalues should approach certain specific functions when the size of matri- 
ces grows. The choice of these functions should depend on the overall symmetry properties of 
matrices such as whether their are complex Hermitian or real symmetric, and on the properties 
of the global limit density of eigenvalues. Roughly speaking, if the global density is positive 
and smooth at a point x, then the limit local behavior at x should be described by the deter- 
minantal point process with the sine kernel. If the global density has a square root singularity 
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at x, then the local behavior around this point should be described by the determinantal point 
process with the Airy kernel, and so on. 

The universality hypothesis originated in pioneering works of Dyson and Mehta. (See, for 
example, the discussion in [16].) It was verified in various particular cases in ||T9l , ||20l , 
10. Very recently, universality was proved for the important cases of Wigner matrices and 
beta ensembles. See Q, l22l and a review in JSJ. 

Now suppose that A and B are different from a multiple of the identity matrix. Is it true 
that the behavior of eigenvalues of the matrix A + UBU* is universal? 

The goal of this paper is to investigate this question in the case when A and B are random 
orthogonal projections, which we will denote P and Q. That is, P and Q are Hermitian (or 
real symmetric) matrices such that P 2 = P and Q 2 = Q. We will assume that the matrices 
P and Q are N-by-N, and their ranks are p and q respectively. We will always consider 
orthogonal projections in this paper so we omit the adjective orthogonal. 

We will call a projection matrix P a random Hermitian (or real symmetric) projection if 
its entries are complex (or real) random variables and their joint distribution is invariant with 
respect to unitary (or real orthogonal) conjugations. 

In a summary, our main finding is that in most situations the behavior of eigenvalues of 
P + Q for large N is universal in the sense that it is similar to the behavior of the eigenvalues 
of classical matrix ensembles. The exceptional case is the local behavior of eigenvalues of 
P + Q near x = 1 when the ranks of P and Q are both approximately N/2, and we discuss 
this case below. 

The reason why the eigenvalue distribution of P + Q is easier to analyze than that of the 
general case A + B is that there is a relation of this distribution to the distribution of the Jacobi 
ensemble of random points. We formulate this relation in a slightly greater generality. 



Theorem 1.1. Suppose P is a random Hermitian (or real symmetric) projection of rank p, 
and Q is an independent Hermitian (or real symmetric) random projection of rank q, p < q. 
Assume p + q < N and 8^0. Then with probability 1, P + 6Q has a = q — p eigenvalues 
8 and b = N — p — q eigenvalues 0. The density of the absolutely continuous part of the 
eigenvalue distribution of P + 6Q is induced by a transformation of the following densities: 
(i) 

p 

fW(x u ...,x p ) = c 2 Yl \xi — Xj\ 2 JJ (1 — Xi) b , Xi G [0, 1] , (1) 

l<i<j<p 8=1 

in the Hermitian case, and 
(ii) 

/«(xi,...,x p ) = ci n \xi-xAjJxf (1- Xi ) b ' , Xi E [0,1]. (2) 

l<i<j<p 1=1 
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in the real symmetric case, with a' = (a — 1) /2 and b' = (b — 1) /2. 
The transformation is two-fold and given by the formula 

\ = ~(i + e ± ^(i-#) 2 + W) . 

This theorem allows us to reduce the study of spectral properties of P + Q to the study of 
the corresponding Jacobi ensembles. In particular, consider a sequence of matrices Pn, Qn 
with increasing size N. We use notation pn, qN to denote their ranks. 

Corollary 1.2. LetAf(I) denote the number of eigenvalues of the Hermitian matrix Pn + Qn 
in an interval I. Assume that pn /N — > p > and qN/N^-q>pasN^-oo, and suppose 
that p + q < 1. If I C (0, 1) or I C (1, 2) arcc? I has positive length, then Var (Af (I)) ~ 
7r~ 2 log N for large N. Moreover, the normalized random variable 

Af (/) - EAf (I) 
VVar(AT(I)) 

converges to a standard Gaussian random variable. 

Note: The symbol C denotes proper inclusion. A similar result holds for real symmetric 
matrices Pn + Qn with a different value of the limit variance. 

Proof of Corollary: Let for definiteness J = (a, j3) C (0, 1) and / = (1 + a 2 , 1 + /3 2 ) C 
(1, 2) . Then M (I) = M' (J), where M' {J) denote the number of points from ensemble CO 
in interval J. This ensemble is determinantal, hence we can apply the theorem due to Costin, 
Lebowitz, and Soshnikov ((H and 11211 ) and conclude that 

N' (J) - W (J) 
VVar (AT' (J)) 

converges to a standard Gaussian random variable, provided that Var (TV"' (J)) grows to infin- 
ity. 

Under the assumption that pn/N — > p > and qN /N — > q > pas TV — > oo, and that 
p + q < 1, it is known that the limit kernel for the corresponding Jacobi ensemble (Q]) is 
the sine kernel. See, for example, [j4] or [Q]. Hence the calculation of the variance can be 
done similar to the calculation of variance in the case of the Gaussian Unitary Ensemble (see 
calculations in ifTTl and related formula (16. 1.3) and Appendix A.38 in [ 16]). This calculation 
results in the variance asymptotically equal to it~ 2 log N for large N. □ 

Similarly, the other properties of eigenvalues of the Hermitian matrix Pn + Qn in the 
interval J are analogous to the properties of the points from the Jacobi ensemble for the simple 
reason that the transformation A— > (X — l) 2 is locally linear on this interval. In particular, in 
the limit these properties can be described by the sine kernel. 

For the behavior at the edge of the support, the results are also similar to the case of the 
classical ensembles and this behavior can be described by either the Bessel kernel (if p+q = 1) 
or the Airy kernel (if p + q < 1). The theorem below gives details for the case p + q < 1. 
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Theorem 1.3. Let Pn and Qn be two random projections in C with ranks pn and qN 
respectively. Assume that Iiuin^oo = P and 'limjv^oo = q, and suppose thatp+q < 1. 



where F2 (u) is the Tracy-Widom distribution function for the edge of the Gaussian unitary 
ensemble. 

This theorem directly follows from recent results about the distribution of the largest eigen- 
value in the Jacobi ensemble with changing parameters. A similar result holds for real sym- 
metric projections except that a is different and the probability converges to F\ (u) instead of 



The only unusual situation occurs at x = 1 when this point belongs to the support of the 
absolutely continuous part of the limit eigenvalue distribution, which occurs when p = q = 
1/2. We will analyze this special situation in the following theorem. 

Let J a (z) denote the Bessel function with index a. 

Theorem 1.4. Let Pn and Qn be two random projections in C N with ranks pn and qN 
respectively. Assume that qN — PN = a > 0, and N — pn — qN = b > 0. Let sn = 
t/ (\/2pAf) , where t > 0. Then the expected number of eigenvalues of P + Q in the interval 
I N = [1,1 + s N ] is 




li = y/q(l-p) + y/p{l - q). 



Then there exists a > such that for every u 




F2 (u) . 



EM(I N )= [ x ( J a {xf - J a+1 (x) J a _! (x)) dx + o(l) 



Jo 



for large N. 



By using the power series for the Bessel functions, 




we obtain the following corollary. 



Corollary 1.5. For large N, 




where R (t) is a dijferentiable function such that lim^o R(t)/t 2a+A exists and finite. 



(Here the asymptotic equivalence ~ is understood relative to increase in N.) For example, 
if a = 0, then EM (I N ) ~ t 2 /2 + R(t) where lim^ R (*) A 4 < 00. 
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In order to appreciate the significance of this corollary, note that under assumptions of 
Theorem 11.41 the eigenvalue distribution of Pn + Qn weakly converges (in probability) to a 
fixed probability distribution which is absolutely continuous on (0, 2) and has the density 



n \A( 2 — ^) ' 

This can be shown either by using a relation of this ensemble with free probability theory (see 
lfl8l ) or by computing the limit distribution for the corresponding Jacobi ensemble. 

In particular, the limiting density is not zero at A = and does not depend on a or b. Hence 
we can see that the expected number of eigenvalues in the interval In cannot be estimated as 
N x length of In x density + terms of smaller order, which would result in a term of the order 
ct with c > 0. For the ensemble Pn + Qn, the expected number of eigenvalues in In is about 
t 2 /2 and significantly smaller than this estimate. 

Thus, the local behavior near x = 1 differs strongly from the local behavior of the GUE 
eigenvalues. In the case of GUE, a similarly scaled limit gives KAf (In) ~ p{x)t where p (x) 
is the density of the semicircle law. Hence, the result in Corollary 1 1.5 l is in contradiction with 
the heuristic that the local behavior of eigenvalues of an Hermitian ensemble is similar to that 
of GUE eigenvalues provided that the limit eigenvalue density is non-zero and smooth at a 
given point. (For example, at the beginning of Section 6.3.3 on p. 113 of 03), it is stated that 
"In the context of Hermitian matrix models, the universality may be stated as the fact that the 
global eigenvalue regime determines the local eigenvalue regime.") 

A possible intuitive explanation of this phenomenon is the additional symmetry of the 
model. This symmetry forces eigenvalues either land exactly on i = 1, or come in pairs 
located symmetrically around x = 1 . In both cases the repulsion between eigenvalues ensures 
that for small t the interval (0, i/2pjv) has fewer eigenvalues than the density function predicts. 

Another interesting feature of the result in Theorem ll.4l is that the limit does not depend 
on b, that is, on the excess of the dimension N of the ambient space over the total rank of P 
and Q. The only provision is that b stays constant as N grows. 

Hence the local behavior of eigenvalues near x = 1 exhibits universality with respect to 
parameter b but not with respect to parameter a. The difference between ranks of P and Q 
does matter for the limiting local behavior of eigenvalues near x = 1 despite the fact that it 
does not influence the global limiting distribution. 

For the case of real symmetric matrices we have an analogous result. 

Theorem 1.6. Let Pn and Qn be two random projections in M, N with ranks pn and qN 
respectively. Assume that qN~PN = a> 0,and N—pN—qN = b> O.LetSN = t/ (\/2pAr) j 
where t > 0. Then for large N the expected number of eigenvalues of P + Q in the interval 
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I N = [1,1 + s N ] is 

E7V (I N ) = X (^Ja'+l (xf - J a '+2 (x) J a' 0)) dx 





+ f J a '+l(x)[j J„<_j (//)(///- 1 ) d.r + o{l). 



where a' = (a — 1) /2. 



By using the power series expansion for the Bessel functions we obtain the following 
consequence. 

Corollary 1.7. For large N, 

EAA (I N ) ~ 1 +R( t ). 

where R (t) is a differentiable function such that lim^o R{t) /t^ 7 ^ 2 exists and finite. 
For example for a = 0, MjV (I N ) ~ 1/ (3^/^/2) t 3 / 2 + P(i), where lim t _> A ?/2 < 



00. 

The rest of the paper organized as follows: we prove Theorem ll.ll in Section[3]after proving 
some results about the relation of eigenvalues of P + Q and of PQP in Section [2] Section [4] 
is about the (universal) edge behavior of eigenvalues. Section [5]proves Theorems 11.41 and 1 1 . 6 1 
about the exceptional case when the behavior of eigenvalues is not universal. And Section [6] 
concludes. 

2. Relation of eigenvalues of P+Q and PQP 

Lemma 2.1. Let P and Q be orthogonal projections and assume that RangePnRangeQ = 0. 
Let ti > be an eigenvalue of operator PQP. Then, 



\ = \ f i + e±\/(i-e f -!- wt, 



are eigenvalues of P + 



Proof: Let n, be an eigenvector of PQP with eigenvalue ti and define Vi := QPu{. Then, 
we have 

QPui = Vi, and Pvi = Uui. 

Clearly, Uj £ RangeP, Vi 6 RangeQ. Hence, Pu.j = Ui and Qvi = Vi. 
Define 



<* = £~ ( - 1 + ± V (l " e ? + m i ) ■ 



1 

2t 
Then, 

(P + 9Q) (ui + c^u;) = Ui + aiUm + Ovi + 9a.iVi = (1 + ( m + —77-7—^°*^ 

V Otiil + Ctiti) 

= \i (Ui + OliVi) . 
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Note that Uj + otiVi ^ by the assumption that RangeP n RangeQ = 0. Hence A, is an 
eigenvalue of (P + 6Q).U 

Lemma 2.2. Suppose P is a random projection of rank p, and Q is a projection of rank q, 
p < q. Assume p + q < N and 9^0. Then, with probability 1, the matrix P + 9Q has 2p 
distinct eigenvalues 

A± = Ul + 0±^(l-0) 2 + W) 

where ti,,.,,t p are non-zero eigenvalues of operator PQP, ti £ (0, 1) . In addition, with 
probability 1, P has q — p eigenvalues 9 and N — p — q eigenvalues 0. 

Proof: Let [L] denote the dimension of linear space L. Then, by using the fact that with 
probability 1 , the range and kernel spaces of P and Q are in general position relative to each 



other, we can calculate: 

[KerPnKerQ] = N - p - q, (3) 

[RangeP n RangeQ] = (jp + q- N) + = 0, (4) 

[KerP n RangeQ] = (N-p + q- N) + = q-p, (5) 

[RangeP n KerQ] = (p + N — q — N) + = 0, (6) 

[KerQP] = [RangeP n KerQ] + [KerP] (7) 

= N-p, (8) 

[RangeQP] = N - [KerQP] = p, (9) 



which implies that QP has p positive singular values, and therefore PQP has p positive 
eigenvalues. In addition, (01) implies that all singular values of QP are smaller than 1. 

Therefore, by previous lemma Q + P has at least 2p eigenvalues (possibly counting with 
multiplicities) that are distinct from and 9. 

With probability 1 the singular values of QP in the interval (0, 1) are distinct. Indeed, the 
set X of P such that PQP has a multiple eigenvalue t with t G (0, 1) is characterised by an 
algebraic condition and therefore it is algebraically closed. The set of all P is algebraically 
irreducible. Hence X is either coincide with the set of all P or has a smaller dimension than 
the set of all P and therefore has measure zero. It is possible to construct an example such 
that all eigenvalues in the interval (0, 1) are distinct. Therefore the measure of X is zero. 

Therefore, the corresponding eigenvalues of Q + P are distinct with probability 1. 

In addition [KerP n KerQ] is an eigenspace of Q + P with dimension N — p — q, and 
eigenvalue 0, and [KerP n RangeQ] is an eigenspace of Q + P with dimension q — p and 
eigenvalue 9. By counting multiplicities we conclude that we found all eigenvalues. □ 

3. Relation of eigenvalues of P+Q and the Jacobi ensemble 

By using known results about the eigenvalues of PQP, we are able to prove Theorem ll.il 
(These eigenvalues have been studied in [4] in connection with free probability.) 
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Proof of Theorem ll.lt By a change of coordinates we can assume without loss of gener- 
ality that Q is a diagonal matrix that has Q„ = 1 for i < q and Q„ = for i > q. 

We will further assume that P is a projection on the span of p independent random vectors 
v\,...,v p . The components of a vector Vi are independent standard Gaussian variables (real 
in the symmetric case and complex in the Hermitian case). If X is an N-by-p matrix with 
columns v%, . . . , t> p ,then we can write: 



X 



A 
B 



where A is a q-by-p matrix and B is a (N — q)-by-p matrix. After orthonormalization, we get 

-1/2 



Y 



and the projection on the span of v\ , 



A 
B 



AA + B'B 

,v p is 
P = YY'. 



(10) 



Formula (fTOl ) implies that the non-zero q-by-q block of QPQ equals A (a'A + B'B) A'. 
Hence, the non-zero eigenvalues of PQP correspond to non-zero eigenvalues of a p-by-p 
matrix M = ^A' A + B'B^j A'A It is known that the joint distribution of eigenvalues 
of M follows the law of the so-called Jacobi ensemble with parameters (p, N — q, q). (See 
Section 3.6 and 7.2.5 of iflOl for more information about the Jacobi ensemble.) Namely, the 
density for the complex case is 

p 



xi e [0, l] 



/ (2) (Xl, . . . ,X p ) = C 2 ]"[ \Xi - Xj\ 2 Y\_X°; 2 (1 - Xi 

l<i<j<p i=l 

where a>2 = q — p and 62 = N — p — q (see Section 8 in [ 12]). 
The density for the real case is 

1 

fW{x 1 ,...,x p ) = c 1 H \xi - Xj\Y[x^ (1 - Xi ) bl , Xi G [0, 1] 

l<i<j<p i=l 



where a\ = (q — p — 1) /2 and b\ 
completes the proof. □ 



(JV - q - p - 1) /2, . (see Muirhead, Thm 3.3.4). This 



4. Eigenvalues at the edge 

In the first step of this section we establish the limit of the eigenvalue distribution on the 
global scale. In particular we identify the support of the limit distribution. Then we study 
the behavior of the eigenvalues at the edge by using the connection between eigenvalues of 
Pn + Qn and the Jacobi ensemble. We rely on known results about the edge behavior of 
eigenvalues from the Jacobi ensemble. As an aside, we derive the limit global density for the 
Jacobi ensemble with variable parameters. 
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We abuse notation a bit and denote lim.jv->oo by p and limjv->oo ^ by q. That is, in 
this section p and q denote the limits of fractional ranks of matrices Pn and Q n ■ We also use 
the notation 

a : = q -p= hm — , 

N-^oo N 

b : = 1 — p — q = lim — , and 

A-^oo N 

c : = p(l-q) + q(l-p) = (l + a 2 -b 2 )/2. 

Let Aj be ordered eigenvalues of Pjy + Qn, counted with multiplicity. Let 

JV 



MAT- 

t=l 

that is, /iiv is the eigenvalue distribution of P/v + Qn- 

Proposition 4.1. WiY/z probability 1 f/je eigenvalue distribution of Pn + Qn weakly converges 
to a probability measure p. The absolutely continuous part of this measure is supported on 
I\ U I2, where 11 = 1 + 1,12 = 1 — I, and 



I 



Vg(i-p)- + 



r/je density of this part is given by 



1 A/-a 2 + 2c(x-l) 2 -(x-l) 4 

H{ ' TT x(x-i)(x-2) 

In addition, p has an atom with weight \a\ at x = l.Ifb > then fx has an atom with weight 
b at x = 0; and ifb<0 then p has an atom with weight —b at x = 2. 

Proof: This result is a particular case of a theorem in [18 1, with the limit distribution 
computed by using the standard techniques from free probability. □ 

As a consequence, we can write down the limit level distribution for the Jacobi ensemble 
with a changing weight function. (This limit level distribution is known and goes back to the 
research of Wachter in [23 1 . However, we obtain it by using free probability calculations and 
this seems to be knew.) 

Let 

N 

f®( Xl ,...,x N ) = c 2 J] \ Xi - Xj \ 2 J{xl N {l- Xl ) tN , Xi e [0,1], (11) 

^<i<j<N i=l 

Recall that the level density of this distribution is simply the marginal distribution: 

p$ (x; s N ,t N ) = / / (2) (x, x 2 , . . • , x N ) dx 2 ... dx N . 
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(2) 

Proposition 4.2. Let p N [x; sjv,ijv) is the level density of the Hermitian Jacobi ensemble 
with parameters sn and i/v- Assume that sn/N — > s > and t]y/N — > t > as N — > oo. 
Define: 



2 + t + s' 
1 ^ g 2 -t 2 

(2 + t + s) 



- ( 1 + — 2 I ' an d 



d ■■ =7^- W(l + *)(! + *) (! + * + «)■ 

(2 + t + s) 

r/zen t/ze feveZ density pj^ (x; sn^n) converges to the limit p^ (x; s, t) which is supported 
on the interval 

I = [c — d, c + d] . 

Its density is 

(2) , 2 + t + s \/-a 2 + 2cx - x 2 
9 {X} = 2^ x(x-l) ■ 

Proof: We use the correspondence between the Jacobi ensemble and the eigenvalues of the 
sum of random projections + Qn established in Theorem ll.il For this we use projectios 
P/v and Qn with ranks p^ and and make sure that pn/N and /N approach p and q, 
which are defined as follows: 

1 , 1 + 8 

p := and q :- 



2+t+s 2+t+s 
By using Theorem 11.11 it is easy to check that the limit eigenvalue density of this matrix en- 
semble must coincide with the limit level density of the Jacobi ensemble. Then, the conclusion 
of the theorem follows from Proposition 14.11 □ 

Now let us find out what is the edge behavior of eigenvalues of P/y + Qn- The edge 
behavior of the Jacobi ensemble with variable parameters was analysed in |4l and lTT4l . (The 
case with the fixed parameters was studied in [9 ] and |[T71 .') 

In particular, Johnstone found that the distribution of the largest eigenvalue is concentrated 
near a point x, which he defines as follows. 
Let 

t—s , t+s 
cos <p = and cos 7 



2+t+s 2+t+s 
Then, 

1 — COS ((ft + 7) 
X = 2 • 

(We change Johnstone's notations a bit to adjust them to our situation since we consider the 
Jacobi ensemble on [0, 1] instead of [—1,1]. In addition, we omitted some terms of order 
O (N -1 ) which are important for statistical applications but irrelevant from the asymptotic 
point of view.) 

An easy verification shows that x coincides with the upper bound of the limiting distribu- 
tion support which we derived in Proposition 14.21 that is, x = c+ d. 
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Johnstone proved the following result (see Theorem 3 in lfl4l and a related result in Theo- 
rem 4.15 of Hi): 

Theorem 4.3 (Johnstone). Assume that t > and let xn\ denote the largest point from N 
points sampled according to distribution Alii . Let sn/N — > s > and t^/N — > t > as 
N —7- oo. Then, there exists a > 0, such that for every u, 

[ tanh -1 (xm) — tanh -1 (x) ) 

p { 4iF575 £ "} ^ <"> • 

where F<i (u) is the Tracy-Widom distribution function for the edge of the Gaussian unitary 
ensemble. 

Johnstone gives an explicit expression for a and an explicit bound for the error in this 
convergence result. He uses a slightly different definition of x which improves the second 
order convergence properties. Also note that while the inverse hyperbolic tangent improves 
the convergence properties for relatively small N, it is irrelevant from the asymptotic point of 
view. A similar theorem was proved by Johnstone for the real symmetric Jacobi ensemble, 
with a different a and F\ (u) instead of F2 (u) . 

Proof of Theorem 11.3 1 This result follows immediately from Theorems 11.11 and |4. 3 1 □ 

5. Eigenvalues of P+Q near x=1 
5.1. Hermitian case. The density 

n 

f(x 1 ,...,x n ) = c ] j \xi - xj\ 2 Y\_Xi (1 - Xi) b , Xi e [0,1] (12) 

l<i<j<n i=l 

can be written as a determinant: 

f(x 1 ,...,x n ) = j^det(K n (xi,Xj)) i)j=1 _ n ; 
n-l 

K n (x, y) = (x) Q k (y) 

fc=0 

(V) - Qn-l (x) Q n (y) 

— C-n 1 

x-y 

where Qk (x) are orthonormal polynomials with respect to the weight x a (1 — x) b on the 
interval [0, 1] and c n is an appropriate constant. 



Lemma 5.1. We have 



lim -^K n ( 



>oo2n 2 \2n 2, 2n 2 J 2 
where J a (t) denote the Bessel function of index a 



- (j a (tf - J a+1 (t) J a _i (t)) , 
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Proof: This is a consequence of a remark in Section 7.2.5 and Proposition 7.2.1 in iflOl . □ 
Proof of Theorem ll.4t By using the relation between the eigenvalues of Pjy + Qn and 
of the Jacobi ensemble, we have EA/" (In) = KM ; (Jn) , where KM' (Jn) is the expected 
number of eigenvalues of the Jacobi ensemble in the interval Jn = [0, tQ/(2p N )]. Hence, by 
using Lemma IBTTl we can write: 

/■*§/(2pSr) 

EM (In) = / K PN (x,x)dx 
Jo 

*° * ( _e_ t 2 \ 



to 



t (Ja (tf ~ Ja+1 (t) Ja-1 (t)) dt + o(l). 



□ 



5.2. Real symmetric case. In the real symmetric case the density of eigenvalue distribution 
and all correlation functions (i.e., the marginals of this density) can be written using quaternion 
determinants: 

f(x 1 ,...,x n ) = —Qdet(K n (x i ,x j )) ij=ln , 
where K n (xj , xj ) is a quaternion function. The matrix representation of this quaternion is 

f S n (x,y) I n (x,y) \ 
{ D n (x,y) S n (y,x) ) ' 



K n (x, y) 
where 



d 

D n (x,y) = —S n (x,y), 

f v 1 

I n (x,y) = -J S n (x, z) dz + -sgn (y - x) , 

and S n (x,y) can be defined in terms of skew-orthogonal or orthogonal polynomials. (See 
Propositions 6.3.2, 6.3.3, and 6.4.3 in IflOl .) 

The expected number of eigenvalues in an interval / for the Jacobi ensemble can be com- 
puted as 

J Sp N (x, x) dx. 

The scaling limit for S n (x, y) at the edge of the spectrum is given in Section 7.9.1 of IflOl . 
We change the formula a bit since we use interval [0, 1] instead of [—1, 1] 



'Ml 

where 



/ ,/„ +l . M r/r). 



K hard (2 2\ _ xJ a+l ( x ) J a (y) ~ V J a+l (v) J a ( x ) 

\*>V)- (x 2 -y 2 ) 
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In particular, 

S hard (x\ X 2 ) = I (j a+1 (^) 2 - J a+2 (x) J a (X)\ (l " [ Ja + 1 (V) dv 



Hence, we can write 



¥J\f (In) = I S PN (x, x) dx 



o 



o »^?3j\r I „ o 1 „ o I Oil 



2 "PJV I 9 2 ' o„2 

o Pn \ z Pn z Pn 

to 



t {Ja+l {tf ~ Ja+2 (t) Ja (t)) 
+ ° J a+1 (t) (l - J J a+1 (v) dv} dt + o(l). 



□ 



6. Conclusion 



We showed that the local behavior of eigenvalues of the sum of two random orthogonal 
projections satisfy the universality paradigm in general, with a single exception. The excep- 
tion can perhaps be explained by an additional symmetry present in the model. 

This evidence give some support to the conjecture that the local behavior of eigenvalues of 
the sum of two non-scalar large random matrices is universal. 
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